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ABSTRACT 

We  develop  parametric  modeling  and  estimation 
methods  for  STAP  data  based  on  the  results  of 
the  2-D  Wold-like  decomposition.  We  show  that 
the  same  parametric  model  that  results  from  the 
2-D  Wold-like  orthogonal  decomposition  naturally 
arises  as  the  physical  model  in  the  problem  of  space- 
time  processing  of  airborne  radar  data.  We  ex¬ 
ploit  this  correspondence  to  derive  computation¬ 
ally  efficient  parametric  fully  adaptive  and  partially 
adaptive  detection  algorithms.  Having  estimated 
the  parametric  models  of  the  noise  and  interference 
components  of  the  field,  the  estimated  parameters 
are  substituted  into  the  parametric  expression  of 
the  covariance  matrix  to  obtain  an  estimate  of  the 
interference-plus-noise  covariance  matrix.  Hence 
the  fully-adaptive  weight  vector  is  obtained.  More¬ 
over,  it  is  proved  that  it  is  sufficient  to  estimate 
only  the  spectral  support  parameters  of  each  inter¬ 
ference  component  in  order  to  obtain  a  projection 
matrix  onto  the  subspace  orthogonal  to  the  inter¬ 
ference  subspace.  The  proposed  partially  adaptive 
parametric  processing  algorithm  employs  this  prop¬ 
erty.  The  proposed  parametric  interference  mitiga¬ 
tion  procedures  can  be  applied  even  when  only  the 
information  in  a  single  range  gate  is  available,  thus 
achieving  high  performance  gain  when  the  data  in 
the  different  range  gates  cannot  be  assumed  sta¬ 
tionary. 

1.  INTRODUCTION 

We  propose  a  new  approach  for  parametric  modeling  and 
estimation  of  space-time  airborne  radar  data,  based  on  the 
2-D  Wold-like  decomposition  of  random  fields.  The  goal  of 
space-time  adaptive  processing  is  to  manipulate  the  avail¬ 
able  data  to  achieve  high  gain  at  the  target  angle  and 
Doppler  and  maximal  mitigation  along  both  the  jamming 
and  clutter  lines.  Because  the  interference  covariance  ma¬ 
trix  is  unknown  a  priori,  it  is  typically  estimated  using  sam¬ 
ple  covariances  obtained  from  averaging  over  a  few  range 
gates.  Next,  a  weight  vector  is  computed  from  the  inverse 
of  the  sample  covariance  matrix,  [l]-[5].  In  [8],  an  approach 
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that  bypasses  the  need  to  estimate  the  covarieince  matrix 
was  presented:  The  data  collected  in  a  single  range  gate 
was  employed  to  obtain  a  least  squares  estimate  of  the  sig¬ 
nal  power  at  each  hypothesized  DOA,  through  evaluation 
of  a  weight  vector  constrained  to  null  the  unknown  inter¬ 
ference  and  noise.  In  [9]  a  simple  ad-hoc  model  of  the  clut¬ 
ter  signal  and  covariance  matrix  is  proposed.  The  model 
represents  the  spectral  density  of  the  clutter  as  a  sum  of 
Gaussian-shaped  humps  along  the  support  of  the  clutter 
ridge.  In  [10]  this  model  is  employed  to  estimate  the  clut¬ 
ter  covariance  matrix  from  the  data  observed  in  a  single 
range  gate. 

In  this  paper,  we  suggest  to  adopt  the  2-D  Wold-like  de¬ 
composition  of  random  fields,  [6],  as  the  parametric  model 
of  the  observed  data.  Employing  this  model,  we  derive 
computationally  efficient  algorithms  useful  for  parametri¬ 
cally  estimating  both  the  jamming  and  clutter  fields.  The 
estimation  procedure  we  propose  is  capable  of  producing  es¬ 
timates  of  the  interference  signals  parametric  models  even 
from  the  information  in  a  single  range  gate.  Hence,  no  av¬ 
eraging  over  a  few  range  gates  is  required,  achieving  high 
performance  gain  in  the  practical  case  when  the  data  in  the 
different  range  gates  is  non-stationary.  Having  estimated 
the  interference  terms  parametric  models,  their  covariance 
matrix  can  be  evaluated  based  on  the  estimated  parame¬ 
ters.  Moreover,  the  problem  of  evaluating  the  rank  of  the 
low-rank  covariemce  matrix  of  the  interference  is  solved  as  a 
by-product  of  obtaining  the  parametric  estimates  of  the  in¬ 
terference  components-  Once  the  parametric  models  of  the 
interference  components  have  been  estimated,  several  alter¬ 
native  detection  procedures  are  available.  In  this  paper  we 
present  two  such  methods:  the  parametric  fully-adaptive 
processing,  and  the  parametric  partially-adaptive  process¬ 
ing. 

2.  THE  RANDOM  FIELD  MODEL 

In  this  section  we  shall  briefly  describe  the  2-D  Wold-like 
decomposition  of  random  fields,  [6].  Let  {y(n,  m)},  (n,  m)  € 
Z^y  be  a  complex  valued,  regular,  homogeneous  random 
field.  Then,  y(n,  m)  can  be  uniquely  represented  by  the 
orthogonal  decomposition 

2/(n,  m)  =  iy(n,  m)  -h  v(n,  m)  .  (1) 

The  field  {a;(n,m)}  is  a  deterministic  random  field.  The 
field  {w{ny  m)}  is  purely-indeterministic  and  has  a  unique 
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Figure  1:  RNSHP  support;  example 
with  a  =  2  and  =  1. 


white  innovations  driven  moving  average  representation, 
given  by 

w{n,  m)  =  ^  b(k,  i)u{n  -kyTn  —  i)  (2) 

(0,0):<(fc.O 

where  S(o,o):<(*.<)  <  0°:  K0,0)  =  1,  and  {ti(n,m)} 

is  the  innovations  field  of  {y(n,  m)}.  The  notation  ^  implies 
that  the  summation  is  performed  over  all  the  samples  that 
are  in  the  of  the  (n,  m)  sample,  where  the  past  is 

defined  with  resF>ect  to  any  selected  choice  of  NSHP  total- 
ordering  on  the  2-D  lattice.  (See,  for  example,  Fig.  1.) 
It  is  possible  to  define.  [6],  a  feimily  of  NSHP  total-order 
definitions  such  that  the  boundary  line  of  the  NSHP  has  a 
rational  slope.  Let  o  and  8  be  two  coprime  integers,  such 
that  a  7^  0.  The  angle  0  of  the  slope  is  given  by  tan  0  = 
/?/a.  (See,  for  example.  Fig.  1.)  A  NSHP  of  this  type  is 
called  rational  non-syrnmetncal  half-plane  (RNSHP).  For 
the  case  where  0  =  0  the  RNSHP  is  uniquely  defined  by 
setting  8=1.  (For  the  case  where  0  =  0  the  RNSHP  is 
uniquely  defined  b>  setting  a  =  1.)  We  denote  by  O  the 
set  of  possible  RNSHP  definitions  on  the  2-D  lattice, 
(Le.,  the  set  of  all  NSHP  definitions  in  which  the  boundary 
line  of  the  NSHP  has  a  rational  slope).  The  introduction 
of  the  family  of  RNSHP  total-ordering  definitions  results  in 
the  following  countabl>‘  infinite  orthogonal  decomposition 
of  the  deterministic  component  of  the  random  field: 

t(n,i7j)  =  p(n,m)  4-  ^  e(a,j9)(n,m)  .  (3) 

(q.5)€0 

The  random  field  {p(n,m)}  is  called  half-plane  determinis¬ 
tic.  The  field  {c(o,d)(«i^)}  is  the  evanescent  component 
that  corresponds  to  the  RNSHP  total-ordering  definition 
(a,/?)€0. 

Hence,  if  {i/(n,  m)}  is  a  2-D  regular  and  homogeneous 
random  field,  then  y{n,  m)  can  be  uniquely  represented  by 
the  orthogonal  decomposition 

y{n,  m)  =  tv(n,  m)  +  p(ti,  m)  -h  ^  e^a,^)  •  (4) 

(a./3)€0 


It  is  finther  shown  in  [6]  that  the  spectral  measures  of  the 
decomposition  components  in  (4)  are  mutually  singular,  A 
model  for  the  evanescent  field  which  corresponds  to  the  RN¬ 
SHP  defined  by  (a,  0)  €  O  is  given  by 

e(a.0)(n,^)  =  («.»”)  (5) 

t=l 

(a,/3) 

=  ^  s\‘^’^\jia-m0)ex.p(j2n—^-^^{n0  +  ina)) 
i=l 

where  the  1-D  purely-indeterministic,  complex  valued  pro¬ 
cesses  -  m0)}  and  {s^“’^^(na  -  m0)},  are  zero- 

mean  and  mutually  orthogoned  for  all  i  7^  j.  Hence,  the 
“spectral  density  function”  of  each  evanescent  field  has  the 
form  of  a  countable  sum  of  1-D  delta  functions  which  are 
supported  on  lines  of  rational  slope  in  the  2-D  spectral  do¬ 
main. 

One  of  the  half-plane-deterministic  field  components, 
which  is  of  prime  importance  in  the  STAP  problem  is  the 
harmonic  random  field 

/i(n,  ^  ( j27r(nu;p  -f  mup)  j  (6) 

where  the  Cp’s  are  mutually  orthogonal  random  variables, 
and  (wp,  i/p)  are  the  spatial  frequencies  of  the  pth  harmonic. 

3.  THE  STAP  MODEL  AND  THE  2-D  WOLD 
DECOMPOSITION 

The  random  field  parametric  model  that  results  from  the 
2-D  Wold-like  orthogonal  decomposition  naturally  arises  as 
the  physical  model  in  the  problem  of  space-time  processing 
of  airborne  radar  data.  In  the  latter  problem  the  target 
signal  is  modeled  as  a  random  amplitude  complex  expo¬ 
nential  where  the  exponential  is  defined  by  a  space-time 
steering  vector  that  has  the  target’s  angle  and  Doppler.  In 
other  words,  in  the  space-time  domain  the  target  model 
is  that  of  a  2-D  harmonic  component  similar  to  (6).  The 
purely-indeterministic  component  of  the  space-time  field  is 
the  stun  of  a  white  noise  field  due  to  the  internally  generated 
receiver  amplifier  no^e,  and  a  colored  noise  field  due  to  the 
sky  noise  contribution.  The  presence  of  a  jammer  results 
in  a  barrage  of  noise  localized  in  angle  and  distributed  over 
all  Doppler  frequencies-  Thus,  in  the  angle-Doppler  domain 
each  jammer  contributes  a  1-D  delta  function  located  at  a 
specific  angle,  and  therefore  parallel  to  the  Doppler  axis. 
In  the  space-time  domain  each  jammer  is  modeled  as  an 
evanescent  component  with  (a,/?)  =  (1,0)  such  that  its  1- 
D  modulating  process  is  a  white  noise  process.  The  ground 
clutter  results  in  an  additional  evaoiescent  component  of 
the  observed  2-D  space-time  field.  The  clutter  echo  from  a 
single  ground  patch  has  a  Doppler  frequency  that  linearly 
depends  on  its  aspect  with  respect  to  the  platform.  Hence, 
clutter  from  all  angles  lies  in  a  “clutter  ridge”,  supported 
on  a  diagonal  line  (that  generally  wraps  around  in  Doppler) 
in  the  angle-Doppler  domain.  A  model  of  the  clutter  field  is 
then  given  by  (5)  with  (a,  0)  such  that  tan  0/a  corresponds 
to  the  slope  of  the  clutter  ridge.  Since  the  rational  numbers 
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axe  dense  in  the  set  of  real  numbers,  an  irrational  slope  of 
the  clutter  ridge  can  be  approximated  arbitrarily  close,  by  a 
rational  one.  Hence  any  clutter  signal  can  be  either  exactly 
modeled,  or  approximated  by  an  evanescent  field. 

We  therefore  conclude  that  the  foregoing  derivation  opens 
the  way  for  new  parametric  solutions  that  can  simplify  and 
improve  existing  methods  of  STAP. 

4.  ESTIMATION  OF  THE  COMPONENTS 

PARAMETERS:  PROBLEM  DEFINITION 

We  next  state  our  assumptions  and  introduce  some  nec¬ 
essary  notations.  Let  {y{n,m)},  (n,m)  €  D  where  D  = 
{(i,  j)|0  <  i  <5  —  1,0  <  j  <T  —  l}be  the  observed 
random  field. 

Assumption  1:  The  purely-indeterministic  component 
{w{n,m)}  is  a  zero  mean  circular  complex  valued  random 
field. 

Assumption  2:  The  number  I  ^ 
evanescent  components  in  the  field,  is  a-priori  known.  This 
assumption  can  be  later  relzixed. 

Assumption  3:  For  each  evanescent  field  the 

modulating  1-D  purely-indeterministic  process  is  a 

zero-mean  circular  complex  valued  process. 

Let  y  =  [y(0, 0), . . . , y(0,  T  -  1), . . . ,  j,(5  -  1,  T  -  1)]^,  and 
let  w,  be  similarly  defined.  Let 

. . . ,  -  1)5), 

sl“’^>(a),  s<“-«(a  -p),...,  -  (T  -  l)/3), 

. . . ,  ((S  -  l)a), . . . ,  -  1)£.  -  (T  -  1)/J)](7) 

be  the  vector  whose  elements  are  the  observed  samples  from 
the  1-D  modulating  process  {5-^’^^}.  Define 

[0,  o,  ...,  (r-l)a, 

A  /?  + a,  /3-^(r-l)a, 

(5-1)5,  (5-l)^-ha,  ...,  (5-l)5-h(T-l)ar(8) 

Given  a  scalar  function  /(v),  we  will  denote  the  matrix, 
or  column  vector,  consisting  of  the  values  of  f{v)  evaluated 
for  all  the  elements  of  v,  where  v  is  a  matrix,  or  a  column 


vector,  by  /(v). 

Using  this  notation,  we  define 

(9) 

Thus,  using  (5), 

we  have  that 

g(a,/3)  ^  ^  ^ 

(10) 

where  ©  denotes  an  element  by  element  product  of  the  vec¬ 
tors. 

Note  that  whenever  na—md  =  ka—ilS  for  some  integers 
n,  m,  k,  i  such  that  0<n,  fc<5-l  and  0<m,^<T-l, 
the  same  element  of  appears  more  than  once  in  the 


vector.  It  can  be  shown,  [7],  that  for  a  rectangular  observed 
field  of  dimensions  5  x  T  the  number  of  distinct  samples 
from  the  random  process  that  are  found  in  the  ob¬ 

served  field  is  Nc  —  (5— l)|a|-f-(T'— 1)|/?1  +  1  — (|a|  — 1)(|/?|  — 
1).  This  is  because  Nc  is  the  number  of  different  “columns” 
one  can  define  on  such  a  rectangular  lattice  for  a  RNSHP 
defined  by  (a,  ff).  We  therefore  define  the  concentrated  ver¬ 
sion^  of  to  be  an  Ac  dimensional  column  vector 
of  non-repeating  samples  of  the  process  Thus  for 

any  (a,  jS)  we  have  that 

(11) 

where  is  rectangular  matrix  of  zeros  and  ones  which 

replicates  rows  of 

Note  however  that  due  to  boundary  effects,  the  vector 
is  not  composed  of  consecutive  samples  from  the  pro¬ 
cess  unless  \a\  <  1  or  |/?|  <  1.  In  other  words,  for 

some  arbitrary  a  and  8  there  are  missing  samples  in 
We  note  that  the  covariance  matrix  which  character¬ 
izes  the  process  is  defined  in  terms  of  the  concen¬ 

trated  version  vector  i.e., 

and  not  in  terms  of  the  covariance  matrix  of  the  vector 

singular  matrix,  given  by  . 

Since  the  evanescent  components  axe  mutually 

orthogonal,  and  since  all  the  evanescent  components  are  or¬ 
thogonal  to  the  purely-indeterministic  component,  we  con¬ 
clude  that  r,  the  covariance  matrix  of  y,  has  the  form 

r  =  rp,+  ^  (12) 

(q,/3)€0 

where  is  the  covariance  matrix  of 

Using  (10)  and  (5)  we  find  that 

0  -  (13) 

5.  PARAMETRIC  ESTIMATION  OF  THE 
INTERFERENCE  COMPONENTS 

In  this  section  we  derive  a  computationally  efficient  algo¬ 
rithm  for  estimating  both  the  jamming  and  clutter  fields, 
based  on  the  above  results.  The  proposed  estimation  algo¬ 
rithm  of  the  spectreJ  support  parameters  of  the  evanescent 
field,  a,  /d  and  is  based  on  the  observation  (see  the 

evanescent  field  model  (5))  that  for  a  fixed  c  =  na  —  m/3 
(i.e.,  along  a  line  on  the  seimpling  grid),  the  samples  of  the 
evanescent  component  are  the  samples  of  a  1-D  constant 
amplitude  harmonic  signal,  whose  frequency  is  The 

algorithm  is  implemented  by  the  following  three-step  pro¬ 
cedure: 

In  the  presence  of  an  evanescent  component,  the  peaks 
of  the  observed  field  periodogram  are  concentrated  along 
a  straight  line,  such  that  its  slope  is  defined  by  the  two 
coprime  integers  a  and  Hence,  several  alternative  ap¬ 
proaches  for  obtaining  an  initial  estimate  of  the  spectral 
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support  p>arameters  of  the  evanescent  comp>onent  can  be  de¬ 
rived  by  taking  the  Radon  or  Hough  transforms,  [12],  of  the 
observ^  field  periodogram.  (The  current  implementation 
employs  the  Hough  transform  for  detecting  straight  lines  in 
2-D  arrays).  However,  due  to  noise  presence,  this  estimate 
may  perturbate.  Since  on  a  finite  dimension  observed  field 
only  a  finite  number  of  possible  (a,  /3)  pairs  may  be  defined, 
the  output  of  the  initial  stage  is  a  set  of  possible  (a,  /?)  pairs 
such  that  the  ratio  £  is  close  to  the  ratio  obtained  for  the 
(a,  /?)  pair  estimated  by  the  Hough  transform. 

For  each  possible  (a,  /?)  pair  we  next  evaluate  the  fre¬ 
quency  parameter  of  the  evanescent  component,  As¬ 

suming  the  considered  (a,  yd)  pair  is  the  correct  one,  we 
know  that  in  the  absence  of  background  noise,  for  a  fixed 
c  =  na  —  mf3  along  a  line  on  the  sampling  grid),  the 
samples  of  the  evanescent  component  are  the  samples  of  a 
1-D  constant  amplitude  harmonic  signal,  whose  frequency  is 
]yy  considering  the  samples  along  such  a  line 
we  obtain  samples  of  a  1-D  constcint  amplitude  harmonic 
signal  whose  frequency  can  be  easily  estimated  using 
any  standard  frequency  estimation  algorithm  the  1-D 
DFT), 

The  test  for  detecting  the  correct  (a,/?)  and  is 

then  based  on  multiplying  the  observed  signal  y(n,  m)  by 

exp(— j27r^»^^^3  {nj3  -h  md)),  for  each  of  the  considered  a, 

and  i/,- triplets,  and  evaluating  the  variance  of  this  signal 
along  a  line  on  the  sampling  grid  such  that  c  —  na  —  m/3, 
Clearly,  the  best  estimate  of  a,  /3  and  is  the  one  that 
results  in  minimal  variance  for  the  1-D  sequence,  as  in  the 
absence  of  noise  the  correct  a,  (3  and  result  in  a  zero 
vziriance. 

Having  estimated  the  spectral  support  parameters  of 
each  evanescent  component,  we  take  the  approach  of  first  es¬ 
timating  a  non-parametric  representation  of  its  1-D  purely- 
indeterministic  modulating  process  {s },  and  only  at  a 
second  stage  we  estimate  the  parametric  models  of  these 
processes.  Hence,  in  the  first  stage  we  estimate  the  partic¬ 
ular  values  which  the  vectors  take  for  the  given  real¬ 

ization,  i.e.,  we  treat  these  as  imknown  constants.  The  es¬ 
timation  procedure  is  implemented  as  follows:  Multiplying 

the  observed  signal  i/(n,m)  by  exp(— (nfS  -{-  md)) 
and  evaluating  the  arithmetic  mean  of  this  signal  along  a 
line  on  the  sampling  grid  such  that  c  —  na  —  m^,  we  have 

1  - 
—  V  y{n,m)exp(-j2iT-:^ — ^(n^  +  ma))  (14) 

iVs  -f 

nfi— m/3=c 

where  Ns  denotes  the  number  of  the  observed  field  samples 
that  satisfy  the  relation  na  —  m/3  =  c.  Once  we  obtained 
the  sequence  of  estimated  samples  from  the  1-D  modulating 
process  the  problem  of  estimating  its  parametric 

model  bwomes  entirely  a  1-D  estimation  problem.  Assum¬ 
ing  the  modulating  process  is  an  AR  process,  and  applying 
to  the  sequence  an  AR  estimation  algorithm  (see,  e.^.,  [13]) 
we  obtain  estimates  of  the  modulating  process  parameters, 
as  well. 


Figure  2:  Spectral  density  of  the  observed  field. 


Finally,  it  is  imF>ortant  to  note  that  we  solve  the  diffi¬ 
cult  problem  of  evaluating  the  rank  of  the  low-rank  covari¬ 
ance  matrix  of  the  interference  as  a  by-product  of  obtaining 
the  par£unetric  estimates  of  the  interference  components: 
Denote  the  number  of  evanescent  components  (interference 
sources)  of  the  field  by  Q-  I*  is  then  shown  in  [11]  that 
the  rank  of  the  interference  covariance  matrix  is  given  by 

rank(r)  =  S  E  Kl  +T  \fik\ -  X)  |a/t| 

fc=i  fc=i  fc=:  fc=i 

the  special  case  where  Q  =  1  and  a  =  1  is  the  well  known 
Brennan  rule,  [3],  of  the  rank  of  the  clutter  covariance  ma¬ 
trix. 

6.  PARAMETRIC  FULLY  ADAPTIVE 
PROCESSING 

Having  estimated  the  parametric  models  of  the  purely  in¬ 
deterministic  and  evanescent  components  of  the  field,  the 
estimated  parameters  can  be  sutetituted  into  (12)-(13)  to 
obtain  an  estimate  of  the  interference-plus-noise  covariance 
matrix  r. 

Let  Vt  denote  the  target  steering  vector,  given  by 

vt  —  b(z37t)  {g)  a(i3t)  .  (15) 

Assuming  a  linear,  uniformly  spaced,  sensor  array  and  a 
uniform  CPI  are  employed  in  our  model,  the  spatial  steering 
vector  a(^)  and  the  temporal  steering  vector  b(t:7)  are  given 

b(n7)  =  , 

respectively.  It  is  well  known  (e.p.,  [3],  p.  57)  that  the 
optimum  space-time  filter  is  given  to  within  a  scale  factor 
by 

w  =  r- V, ,  (16) 

The  test  statistic  z(-cs,d)  is  then  given  by 
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Figure  3:  The  test  statistic  z{zu^  i9). 


2(c7,l?)  =  w"(CT,l?)y  =  V,"(c7,J?)(r  ^)"y.  (17) 

Let  Xf  —  (F“^)^y,  We  thus  have 

z{ts,  ■S)  =  vf  (w,j?)x/  =  b"(j^)  (g  a"(i?)x/  .  (18) 

Reorganizing  the  elements  of  Xf  into  a  T  x  S  matrix  ^ 
where  the  elements  of  the  fcth  row  of  ^  are  x/((^  —  1)5  4- 
1) ...  X/  (fc5),  we  conclude  that  for  a  linear,  uniformly  spaced, 
sensor  array  and  uniform  CPI  the  test  statistic  is  given  by 

T  S 

p=l  9=1 

Thus,  i9)  and  ^(p,  q)  are  a  2-D  DFT  pair,  and  the  test 
is  equivalent  to  finding  the  2~D  frequency  where  the  2-D 
DFT  of  ^(p,  q)  is  maximal. 

To  illustrate  the  operation  of  the  propiosed  solution  we 
resort  to  numerical  evaluation  of  some  specific  examples. 
Consider  a  2-D  observed  random  field  consisting  of  a  sum 
of  a  purely-indeterministic  component  (background  noise), 
a  single  evanescent  (interference)  component,  and  three  har¬ 
monic  components  (targets).  The  purely-indeterministic 
component  is  a  complex  valued  circular  Gaussian  white 
noise  field.  The  evanescent  component  spectral  support 
parameters  are  (a,/?)  =  (1,-2),  =  0.  The  modu¬ 

lating  1-D  purely  indeterministic  process  of  this  evanescent 
component  is  a  first  order  Gaussian  AR  process,  such  that 
its  driving  noise  variance  =  2,  and  c^^’'"^^(l)  = 

—0.5.  There  are  three  targets  which  are  located  at  (0.05, 0), 
(0.15,0.15)  and  (—0.25,0.15),  respectively.  The  observed 
field  dimensions  are  48  x  48. 

Let  us  define  the  experimental  variance  of  each  of  the 
field  components  as  Ew  —  w^w  for  the  purely  indetermin¬ 
istic  component;  Ee  =  for  the  evanescent 

component;  and  Ehi,  =  hf  hjt,  =  1,2,3,  for  each  of  the 
harmonic  components,  where  hjk  is  defined  in  the  same  way 


w  and  are  defined.  In  this  example  we  have  ^  = 

6dB,  while  for  the  three  targets  we  have  =  — 12.8dB, 

=  —  14.5dB,  =  — 15dB.  Due  to  the  strong  in¬ 

terference  component,  the  presence  of  the  three  targets  is 
hard  to  detect  in  the  observed  data  whose  power  spectral 
density  is  depicted  in  Fig.  2.  However  these  targets  are 
easily  detected  by  the  test  statistic  zi'UJ,  'd)  depicted  in  Fig, 
3.  In  Fig.  3,  z{is7,'d)  is  depicted  as  a  function  of  the  two- 
dimensional  frequencies,  le.,  angle  and  Doppler. 

7.  PARAMETRIC  PARTIALLY  ADAPTIVE 
PROCESSING 

Recall  that 

(20) 

Having  estimated  and  using  the  algorithm  in 

Section  5,  the  vector  is  known.  Hence,  demodulating 
we  conclude  using  (10)  that  the  demodulated  vector 
which  we  denote  by  is  given  by 

e(a,/3)  ^  ^  ^  ^21) 

Prom  (11)  we  conclude  that  the  covariance  matrix  of  e^^’^^is 
given  by 


(22) 


In  the  following  it  is  proved  that  since  a  and  jS  are 
already  known,  an  orthogonal  projection  matrix  onto  the 
low-rank  subspace  sp>anned  by  the  evanescent  field  covari¬ 
ance  matrix  can  be  found  without  estimating  the  parcimet- 
ric  model  of  the  evanescent  field  1-D  modulating  process, 
and  hence  without  estimating  Moreover  this  result 

enables  us  to  avoid  the  need  in  both  evaluating  the  field 
covariance  matrix,  and  in  employing  a  computationally  in¬ 
tensive  eigenanalysis  to  the  estimated  covariance  matrix. 

More  specifically,  let  us  construct  the  following  orthog¬ 
onal  projection  matrix 


r(a,/3) 


=  ((Aj“’^^)^A^“’^^)  \aJ“’®>)^(23) 


It  is  easily  verified  (by  substitution)  that  is  an  or¬ 

thogonal  projection  onto  the  range  space  of  f  since  for 
any  ST  dimensional  vector  v 

=  f .  (24) 

Also,  (T^^’^^)2  =  and  (T^^’^^)^  = 

Note  that  since  is  a  sparse  matrix  of  zeros  and 

ones  only,  the  computation  of  is  very  simple. 

The  projection  matrix  onto  the  subspace  orthogonal  to 
the  interference  space  is  therefore  given  by  = 

I  —  Hence  by  projecting  the  demodulated  observed 

data  vector  y  =  y  ©  ((dj'^’^^)^)^  onto  the  subspace  orthog¬ 
onal  to  the  interference  subspace,  a  reduced  dimension  data 

vector  given  by  y  =  y  is  obtained,  such  that 
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Figure  4:  Spectral  density  of  the  field  after  being  pro¬ 
jected  onto  the  subspace  orthogonal  to  the  interference 
subspace. 


the  interference  contribution  to  the  observed  signal  is  mit¬ 
igated.  Remodulating  y  by  evaluating  y  0  followed 

by  sequentially  applying  this  procedure  to  mitigate  each  of 
the  interference  sources,  the  detection  problem  is  reduced 
to  that  of  detecting  a  target  in  the  presence  of  background 
noise  only.  Thus,  in  the  special  case  where  the  background 
noise  is  known  to  be  a  white  noise  field,  the  statistical  test 
is  equivalent  to  finding  the  2-D  frequency  where  the  2-D 
DFT  of  the  processed  data  vector  (organized  back  into  a 
2-D  array)  is  maximal. 

As  an  example  consider  the  same  field  as  in  the  previous 
section.  Due  to  the  strong  interference  component,  the  pres¬ 
ence  of  the  three  targets  is  hard  to  detect  in  the  observed 
data  whose  power  spectral  density  is  depicted  in  Fig.  2. 
However  these  targets  are  easily  detected  in  the  processed 
data  as  illustrated  in  Fig.  4.  This  result  is  obtained  with¬ 
out  estimating  the  parametric  model  of  the  evanescent  field 
1-D  modulating  process,  and  hence  without  estimating  the 
interference-plus-noise  covarieince  matrix.  Since  both  the 
estimation  of  the  interference-plus-noise  covariance  matrix, 
as  well  as  its  analysis  are  saved,  the  proposed  parametric 
partially  adaptive  processing  method  is  robust  and  compu¬ 
tationally  attr€w:tive. 
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